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SECOND CLASS PARTICLES AND CUBE ROOT ASYMPTOTICS 
FOR HAMMERSLEY'S PROCESS 

By Eric Cator and Piet Groeneboom 

Delft University of Technology and Vrije Universiteit Amsterdam 

We show that, for a stationary version of Hammersley's process, 
with Poisson sources on the positive x-axis and Poisson sinks on 
the positive j/-axis, the variance of the length of a longest weakly 
North-East path L{t,t) from (0,0) to {t,t) is equal to 2E{t - X{t))+, 
where X{t) is the location of a second class particle at time t. This 
implies that both E(t — X{t))+ and the variance of L{t,t) are of 
order t^^^ . Proofs are based on the relation between the flux and the 
path of a second class particle, continuing the approach of Cator and 
Groeneboom [Ann. Probab. 33 (2005) 879-903]. 

1. Introduction. In an influential paper Kim and Pollard [8] show that 
in many statistical contexts we are confronted with estimators which con- 
verge at rate n^^^ instead of the usual rate n^/^ and that in this situation 
the limit distribution is nonnormal. They call this phenomenon "cube root 
asymptotics." A prototype of such an estimator is the maximum likelihood 
estimator of a decreasing density, which converges locally at rate n^/^ after 
rescaling to the (almost surely unique) location of the maximum of Brown- 
ian motion minus a parabola. The characterization of this limit distribution 
in terms of Airy functions was given in [6]. 

It has been conjectured that the asymptotics for longest increasing sub- 
sequences, which can be analyzed by studying longest North-East paths of 
Hammersley's process, is related to these cube root phenomena in estimation 
theory and, in particular, that it should be possible to derive the asymp- 
totics along similar lines. However, up till now, the cube root limit theory 
for longest increasing subsequences and longest North-East paths has been 
based on certain analytic relations, involving Toeplitz determinants; see, for 
example, [2] and [3]. 
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Fig. 1. Space-time paths of the Hammersley' s process, with sources and sinks. 



In this paper we will work with Hammersley's process with sources and 
sinks, as defined in [4]. We will give a short description here, based on Figure 
1. We consider the space-time paths of particles that started on the x-axis 
as sources, distributed according to a Poisson distribution with parameter 
A, and we consider the i-axis as a time axis. In the positive quadrant we 
have a Poisson process of what we call a-points (denoted in Figure 1 by x), 
which will have intensity 1, unless otherwise specified. On the t-axis (which 
also sometimes will be called y-axis) we have a Poisson process of what we 
call sinks of intensity 1/A. The three Poisson processes are independent. 

At the time an a-point appears, the particle immediately to the right 
of it jumps to the location of the a-point. At the time a sink appears, the 
leftmost particle disappears. To know the particle configuration at time s, 
we intersect a line at time s with the space-time paths. The counting process 
of the particle configuration at time t is denoted by Lx{-,t), where we start 
counting at the first sink on the t-axis up to (0,t), and continue counting 
on the halfline (0,oo) x {t}, so Lx{x,t) equals the total number of sinks in 
the segment {0} x [0,t] plus the number of crossings of space-time paths of 
the segment [0,x] x {t}. 

The total number of space-time paths in [0, x] x [0, t] is called the flux at 
{x,t). It is in fact equal to Lx{x,t). If A = 1, we will denote Li{x,t) just by 
L{x,t), unless this can cause confusion. The flux Lx{x,t) equals the length of 
a longest weakly (North-East) path from (0,0) to (x,t), where "weakly" 
means that we are allowed to pick up either sources from the x-axis or sinks 
from the 4-axis, before we start picking up a-points. To see this from Figure 
1, trace back a longest weakly NE path from {x,t) to (0,0), and note that 
one will pick up exactly one a-point or one source or one sink from each 
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space-time path. Note that, if < x < y, L\{y,t) — L\{x,t) is the number 
of particles (or crossings of space-time paths) on the segment [x,y\ x {t}. 

A heuristic argument for the cube root behavior of the fluctuation of the 
length of a longest weakly NE path for the stationary Hammersley process 
runs as follows. Suppose, for simplicity, that A = 1. A longest weakly NE 
path with length L{t,t) = Li{t,t) can pick up points from either x- or y- 
axis before starting on a strictly NE path to {t,t). Furthermore, let, for 
-t<z<t, 



(1.1) N{z) 
and 

il.2)Mz) 



number of sources in [0, z] x {0}, if z > 0, 
number of sinks in {0} x [0, \ z\], if z < 0, 



length of longest strictly NE path from (2,0) to {t,t), 
if z > 0, 

length of longest strictly NE path from (0, \z\) to {t,t), 
if z < 0. 



Note that the processes At and N are independent and that 

(1.3) L{t, t) = sup{iV(z) +At{z):-t<z<t}. 

The process z t~^^^{N{zt'^^^) — \z\t'^^^}, \z\ < t^^^, converges in the 
topology of uniform convergence on compacta to two-sided Brownian mo- 
tion, originating from zero. As will be shown below, the expectation of 
t~^/^At{zt'^^^) has an asymptotic upper bound (as t — > 00) of the form 

2t2/3 _ 1^1^1/3 _ 1^2^ 

which is seen by taking expectations and optimizing the choice of A in the 
inequality in Lemma 4.1. This suggests that the distance to zero of the exit 
point where the longest path leaves either x- or y-axis cannot be of larger 
order than t^/'^, since otherwise the Brownian motion cannot cope with the 
downward parabolic drift —jz'^, temporarily assuming that the fluctuation 
of At{z) is of order Op{t^^^). The latter fact we know to be true from the 
analytic approach, not used in our probabilistic approach. On the other 
hand, we will derive this in Section 7; see (7.7). 

The limit behavior of the exit point can be compared to the behavior of the 
location of the maximum of Brownian motion minus a parabola, which plays 
a key role in the asymptotics for the cube root estimation theory, mentioned 
above. The crucial difference, however, is that the exit point is the location 
of the maximum of the sum of two independent processes instead of the 
maximum of just one process. 

We note here that the n^^^ convergence in estimation theory (so slower 
convergence than the usual n^/^-convergence) corresponds to the t^/^ order 
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of the distance to zero of the exit point, which (after a time-reversal ar- 
gument, based on Burke's theorem for Hammersley's process) can be cahed 
"super-diffusive" behavior of a second class particle. The slower convergence 
in estimation theory is caused by the fact that the estimators have an inter- 
pretation in terms of the location of a maximum, just as the exit point for a 
longest weakly NE path in Hammersley's process. 

The key relation which allows us to make the heuristic argument above 
rigorous is Theorem 2.1 of Section 2, which, in combination with (time- 
reversal) results of Section 3, tells us that Var(L(t, t)) = 2EZ(t)_|_, where Z(t) 
is the rightmost point where a longest weakly NE path leaves the rr-axis; see 
(3.2), Section 3. It is shown in Section 4 that this implies EZ(t)_|_ = 0(t^/'^). 
In Section 5 we compare longest strictly NE paths with longest weakly NE 
paths and obtain a bound on the difference between the lengths of these 
paths. This allows us to also obtain a lower bound for EZ(t)_|_ in Section 6. 
Finally, in Section 7 we discuss tightness results for the original Hammersley 
process without sources and sinks, in connection with results of Seppalainen 
[9]. 

Our methods heavily rely on the ideas developed in [4], which concern, in 
particular, the difference in behavior below and above the path of a second 
class particle and Burke's theorem for Hammersley's process, which enables 
us to use time-reversal and reflection. 

2. Variance of the flux and location of a second class particle. We will 
need the concept's second class particle and dual second class particle, which 
also play an important role in later sections. A "normal" second class particle 
is created by putting an extra source at (0, 0) (thus effectively removing the 
first sink), and a dual second class particle is created by putting an extra 
sink at (0,0), thus effectively removing the first source. Define X{t) as the 
location at time t of a second class particle in a stationary Hammersley 
process with source and sink intensity equal to 1 (the symmetric case), and 
X'{t) as the location at time i of a dual second class particle for this case. 

As explained in [4], a "normal" second class particle X{t) jumps to the 
previous position of the ordinary ("first class") particle that exits through 
the first sink at the time of exit, and successively jumps to the previous 
positions of particles directly to the right of it, at times where these particles 
jump to a position to the left of the second class particle. The concept 
of a dual second class particle was also considered in [4], but there it is 
seen as a second class particle for the process "moving from left to right." 
Figure 2 shows the trajectories of a second class particle and a dual second 
class particle. Note that we always have X{t) < X'{t), which is evident from 
Figure 2. 

Now consider a stationary Hammersley process with a-intensity 1, source- 
intensity A and sink-intensity 1/A. Fix x,t > and consider the flux L\{x,t). 
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Denote X\{t) and X'-^{t) as the locations at time t of a second class particle 
and a dual second class particle, respectively. We use the subscript A to in- 
dicate that the distribution of the location of the (dual) second class particle 
depends on A. If A = 1, the subscript is suppressed. We have the following 
result: 



Theorem 2.1. 

Yai{Lx{x,t)) = - Ax + ^ + 2AE(x - Xx{t))_^. 

A 

Remark 2.1. A similar relation between the variance of the flux and 
the location of a second class particle has been proved for totally asymmetric 
simple exclusion processes (TASEP) in [5]. 

Remark 2.2. Note that taking A = \/t/x yields 
Yar{Lxix,t)) = 2XE{x - Xx{t))^. 

Proof of Theorem 2.1. For notational clarity we use the four wind 
directions N, E, S and W to denote the number of crossings of the four 
respective sides of the rectangle [0,x] x [0,t] [so Lx{x,t) = N + W]. Clearly, 
S + E = N + W. We also know from Burke's theorem for Hammersley's 
process (see [4]) that N and E are independent, just like S and W. This 
means that 

YaT{Lx{x,t)) =Yar{W + N) 



{X(t)J) ix,t) (X'{t),t) 
« 



o 



o 
o 
o 



-i^ — & — ^—(f- 



♦ — ^ — * — ll* — 



Fig. 2. Trajectories of {X{t),t) and {X'{t),t). 
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= Var(VF) + Var(iV) + 2 Cov{W, N) 
(2.1) = Var(VF) + Var (iV) + 2 Cov{S + E-N,N) 

= Var(VF) - Var(iV) + 2 Cov(5, A^) 

= ^ - \x + 2Cov{S,N). 
A 

We want to investigate Coy{S, N). It turns out that we can do this by varying 
the source-intensity appropriately. For e > 0, we define a source-intensity of 
X + e. The sinks remain a Poisson process with intensity 1/A. We denote 
expectations with respect to this new source intensity by E^. Define 

a„ = E^(iV|5 = n). 

Note that a„ does, in fact, not depend on e, since we condition on the number 
of sources in [0,x], and the sources outside this interval do not influence N. 
Then 



d_ 

Ye 



e=0 



(^(A + £))\ -.(A+.). 



E 



-e ^ 'a 



■n 



de 

A n n! 

71=0 n=0 

\k{NS)-x¥.{N). 
X 



This shows that 
(2.2) 



Cov(iV, S) = ¥.{NS) - E(iV)E(5) 



^1 



£ = 



where we use that E(5) = Ax. 

We will calculate this derivative in the following manner. Fix, indepen- 
dently, a Poisson process of intensity 1 of a-points in (0, oo)^, a Poisson 
process of sources of intensity A on the x-axis and a process of sinks of in- 
tensity 1/A on the t-axis. Now we add an independent Poisson process of 
intensity e to the process of sources. Define as the number of crossings 
of the North-side [i.e., (0,3;) x {t}] for the process with the added sources. 

Note that if we add an extra source at (z, 0), then increases by 1 if and 
only if Xx{t] z) < x, where Xx{t; z) is the location of a second class particle 
at time t, which started at (^,0). We denote X\[t) = X\{t]{)). This means 
that 

(2.3) E(7V,)=E(7Vo) + e rnHx^{t;.)<.})dz + 0{e''). 

•J 
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d 



Cov(N, S) = \^ 
oe 



E(iV,) 

e=0 

/•X 

J 

(2.4) = A r P(XA(t) <x-z)dz 



= A / Y(x-Xx{t)> z)dz 
Jo 

= \E{x-Xx{t))^. 
Combining this with (2.1) gives 

Var(LA(x,t)) = - Ax + ^ + 2AE(a; - Xx{t))^. □ 

A 

Now consider a stationary Hammersley process with source (and sink) 
intensity 1. We denoted the flux of this process at {x,t) by L{x,t), and the 
location of a (dual) second class particle at time t, which started at (0,0), 
by X{t) [resp. X'{t)]. Note that under the map 

{x,t) ^ (x/A, At), 

a stationary process with source intensity 1 gets transformed into a station- 
ary process with source intensity A (and corresponding sink intensity 1/A). 
This rescaling argument shows 

(2.5) XXx{t)=X{t/\) and XX'x{t) = X'{t/X), 

T) 

where = denotes equality in distribution. 

We would like to bound \ai{L\[t,t)) in terms of Var(L(t,t)) in the case 
where A > 1. Using Theorem 2.1 and (2.5), we get, using the inequality 

{A + B)+<A+ + B+, 

Yai{Lx{t, t)) = -Xt + t/X + 2E(At - t/X + t/X - X(t/A))+ 
< (A - 1/A)t + 2E{t/X - X{t/X))^ 
= (A - l/A)t + Var(L(t/A,VA)). 
If we show the intuitively clear result that, for A > 1, 

(2.6) \ai{L{t/X,t/X)) < Var(L(t,t)), 
we have proved that 

(2.7) YaT{Lx{t,t)) < (A - 1/A)t + Var(L(t,t)). 
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We can show (2.6) by noting that Theorem 2.1 for A = 1 is equivalent to 

(2.8) Var(L(2;, t)) = -x + t + 2 T F{X{t) < z) dz. 

Jo 

Define 

(j){x,t)=Yav{L{x,t)). 

Clearly, (j) is symmetric, since the source and sink intensities are equal, which 
gives reflection symmetry of the process. 

Furthermore, (2.8) shows that (p is a continuously differentiable function, 
with 

di(l){x, t) = -l + 2¥{X{t) < x). 

If we can show that ¥{X(t) <t)> 1/2, we would have proved (2.6), since 
di(j){t,t) > for a symmetric function (/) implies that (j){t,t) is increasing in 
t. 

Since reflecting Hammersley's process in the diagonal preserves the dis- 
tribution, while interchanging the trajectories of X' and X, we know that 

F{X{t) >x)= F{X'{x) <t)< F{X{x) < t). 

Choosing t = x, we see that 

F{X{t) >t)< F{X{t) <t)= F{X{t) < t), 

which shows that F{X{t) <t)> 1/2. As noted before, this proves (2.6). 

3. Connection between second class particles and exit points. As has 

already been noted in the Introduction, we can view the flux Lx{x,t) in 
two ways: it is the number of space-time paths in the square [0,x] x [0,t], 
but is also the length of the longest weakly NE path from (0,0) to {x,t), 
where "weakly NE" means that we are allowed to pick up sources or sinks, 
as well as a-points, as long as we are going North-East. To work with this 
latter representation, which we will mainly use in the symmetric case when 
both the source- and the sink-intensity are 1, we define, for —t<z<t, N{z) 
and At{z) by (1.1) and (1.2). Remember that the processes and N are 
independent and that 

(3.1) L{t, t) = sup{N{z) +At{z):-t<z<t}. 

Another important aspect of this representation is the location at which 
a longest path leaves either the x-axis or the y-axis. Define 

(3.2) Z{t) = sup{z e [-t, t] : N{z) + At{z) = L{t, t)} 
and 



(3.3) 



Z'{t) = mf{z £ [-t, t] : N{z) + At{z) = L{t, t)}. 
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Fig. 3. Relating X and Y . 



We call Z{t) and Z'[t) exit points for a longest path, since there exist longest 
paths that leave the axis on (Z(t),0) [or (0, — or on (Z'(t),0) [or 
(0,— From this definition and using the symmetry of the situation, 
we can see that 

(3.4) Z'{t)<Z{t) and Z{t) = -Z'{t). 

We will need another link between the two representations. We have de- 
fined X{t) and X'{t) as the position at time t of a second class particle, 
respectively dual second class particle, that starts at (0,0). Now define 

t-X{t), iiX{t)<t, 
inf{s>0:X(s) >t} -i, if X(i) > t. 



Y{t) 
and 

Y'{t) 



t-X'{t), ifX'{t)<t, 
mi{s>0:X'{s)>t} -t, ifX'{t)>t. 

Since X'{t) > X{t), we have Y'{t) < Y{t). Figure 3 shows the relation 
between X and Y. 

It also shows two relations which we will be important later: 

a<t =^ {X{t) <a} = {Y{t) >t-a} and 

{X' (t) < a} = {Y' {t) > t - a}, 

(3.5) 

a>t =^ {X{t) >a} = {Y{a) <t-a} and 

{X' (t) > a} = {Y' (a) < t - a}. 

Now consider Figure 4. The left picture shows a realization of the Ham- 
mersley process and two longest weakly NE paths, corresponding to Z{t) 
and Z'[t). The right picture shows the same realization, but now reflected 
in the point (^t,^t). 
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Note that the longest paths become trajectories of a second class and a 
dual second class particle in the reflected process, and that Z{t) corresponds 
to Y{t), while Z'{t) corresponds to Y'{t). Burke's theorem in [4] states that 
the reflected process is also a realization of the stationary Hammersley pro- 
cess, so that we can indeed conclude that 

Z{t) = Y{t) and Z'{t) = Y'{t). 

In particular, this means, using Theorem 2.1 and noting that (t — X{t))+ = 
y(t)+, that 

(3.6) Var(L(i,t)) = 2EZ(t)+. 

4. EZ{t)^ is of order 0(<^/^). We wish to control the exit point Z{t). 
We will do this by considering an auxiliary Hammersley process Lx, coupled 
to the original one, by thickening the sources to a Poisson process of intensity 
A > 1 and thinning the sinks to a Poisson process of intensity 1/A. The 
process At then satisfies the following inequality. 

Lemma 4.1. Let A > 1 and define Lx{x,t) as the flux of Lx at {x,t). 
Then, for < z <t, 

Mz)<Lx{t,t)-Lx{z,0). 

Proof. It is clear that a strictly NE path from (2^,0) to {t,t) is shorter 
than a longest weakly NE path from (z, 0) to (t, t), where this path is allowed 
to either pick up sources of Lx on [z, t] x {0}, or pick up crossings of Lx with 
{z} X [0,t]. However, this longest weakly NE path is equal to the number 
of space-time paths in [z,t] x [0,t] of Lx, which, in turn, is equal to Lx{t,t) 
minus the number of sources on [0,z] x {0}. □ 



z'ity> 



<i> 
<> 

i> 
<i> 
<> 

* — 4 
o 
o 



Z{t) t 



lY'(t) 




Fig. 4. Longest path is distributed as trajectory of a second class particle. 
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We can now show the following theorem. We use the notation a{x) < b{x) 
if there exists a constant M such that, for all parameters x, a{x) < Mb{x). 

Theorem 4.2. LetQ<c< t/E2'(t)+. Then 

P(ZW>cEZ(tU)<^^^(-i + -L). 

Proof. Note that, for any A > 1, 

¥{Z{t) > u} = P{3 z > n : N{z) + At{z) = L{t, t)} 

<F{3z>u: N{z) + Lx{t, t) - Lx{z, 0) > L{t, t)} 

= P{3 z > n : N{z) - Lx{z, 0) > L{t, t) - Lx{t, t)}. 

Since Lx{-,0) is a thickening of L(-,0), we get that Nx-i{z) := Lx{z,0) — 
N{z) is in itself a Poisson process with intensity A — 1. This means that 

P{Z(t) > u} < F{NxMu) < Lxit, t) - Lit, t)}. 

To have a useful bound for all < n < |t, we choose A such that 

ENx-i{u) - E{Lx{t, t) - L{t, t)} = {X-l)u-t(^X + j-2^ 

is maximal. This means that we choose 

Some useful elementary inequalities, that hold for all < n < |t, are 

Am < 2, 

(4.1) ENx^_^iu)-E{LxAt,t) - L{t,t)} > ^/t, 

K - l/An < 2u/t. 

Note that, due to (2.7) and (3.6), 

Var{LA„ (t, t) - L{t, t)} < 2(Var{LA„ {t, t)} + Var{L(t, t)}) 

< 8EZ(t)+ + 2t(A„ - 1/Xu) 

< 8EZ(t)+ + Au. 
Now we can use Chebyshev's inequality: 
P{Z(t) > n} < P{iVA„_i(u) < La„ (t, t) - Lit, t)} 

< P{iVA„_i(n) < ENx^^^iu) - nV(8t)} 

+ F{Lx^{t,t) - L{t,t) > EiVA„_i(n) - uV(8t)} 
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64i2(A„ - l)u 



(4.2) < 



+ P{La„ (t, t) - L{t, t) > E{La„ (t, t) - L{t, t)} + nV(8t)} 
^ t2 64t2(8EZ(t)++4u) 



If t > n > |t, we see that 



P{Z(t)>n}<p|z(t)>^t|<i3 + 



where we use (4.2). This means that (4.2) is true for all Q < u <t. The 
theorem now follows from choosing u = cKZ{t)j^. □ 

With this theorem we can show that EZ(t)+ = 0{t^/^). 

Corollary 4.3. Let Z{t)+ and L{t,t) he defined as in (3.6). Then 

EZ{t)+ Var(L(t,t)) 
limsup — ^ — = limsup — < +oo. 

Proof. Using (3.6), we only have to prove the statement for EZ(t)_|_. 
Suppose there exists a sequence t„, j +00 such that 

lim — = +00. 

n-*co tn I 



tl /II 



Using Theorem 4.2, we see that 

Using dominated convergence [note that t^/ (EZ(t„)_|_)^ is a bounded se- 
quence], this shows that 



POO 

/ P{Z(t„)+>cEZ(i„)+}dc^O, 
Jo 



which would imply the absurd assertion that 

□ 

-r 

As a corollary we get the following: 



Z{tn)+ 1 
\EZ{tn)+] 
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Corollary 4.4. Let c> 1. Then 

P{Z(t)>ct2/3}<^. 

Proof. This is an immediate consequence of Theorem 4.2 and the pre- 
vious corollary. □ 

Remark. This result can be compared to a result on transversal fluc- 
tuations of a longest NE path in [7]. He shows that all longest stricly NE 
paths from (0,0) to {t,t) remain in a strip along the diagonal of width t"' , 
with probabihty tending to 1, as t — > oo, for any 7 > 2/3. To this end, he 
uses the analytic results in [2]. Section 3, in combination with Corollary 4.4, 
shows that the transversal fluctuation at time s < t of a longest weakly NE 
path from (0,0) to {t,t) is of order {t — s)^/^. This is due to the fact that 
any longest weakly NE path lies within the reflected trajectories of a sec- 
ond class particle and a dual second class particle; see Figure 4. Our result 
on weakly NE paths implies the same result for strictly NE paths, as the 
following short argument will show: consider the longest strictly NE path 
that is to the right of all other longest strictly NE paths and suppose, at 
time s <t, this path is to the right of the diagonal. Since the sources and 
sinks are independent of this event, we have that, given this event, there is 
still a probability of at least 1/2 that Z{t) > 0. The corresponding longest 
weakly NE path (so the weakly NE path that is most to the right) cannot be 
to the left of the considered longest strictly NE path (because they cannot 
cross), so the order of the fluctuations to the right of a longest strictly NE 
path cannot be higher than the same order for longest weakly NE paths. For 
fluctuations to the left, a similar argument holds. The remark at the end of 
Section 6 discusses the corresponding lower bound result. 

5. Strictly NE paths and restricted weakly NE paths. To get a lower 
bound on EZ(t)_|_, we need to control the difference between a strictly NE 
path and a weakly NE path in a stationary Hammersley process Li (with 
source intensity 1), where the weakly NE path is only allowed to pick up 
sources m an interval [0,6*^/^] x {0}. To do this, we consider another in- 
dependent Hammersley process Lx on [0,t]^ with source intensity A, sink 
intensity 1/A and o-intensity 1; for this process, the sources, sinks and a- 
points are independent of the corresponding processes for Li. Coupled to 
this process Lx, we consider Lq as the corresponding (nonstationary) Ham- 
mersley process that uses the same a-points, but has no sources or sinks. 

We denote Lo{x, t) as the number of particles (i.e., the number of crossings 
of space-time paths) of the Hammersley process without sources or sinks 
with the segment [0,x] x {t}. Note that, for < z <t, 



(5.1) 



^(0) - At{z) = Lo{t, t) - Lo{t -z,t). 
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This follows from the fact that Lo(t — z,t) equals the length of the longest 
(strictly) NE path from (0,0) to {t - z,t). Also note that {N{z):z£ (0,t)} 
[the number of sources of the process Li in the interval (0, t)] and {Lx(t, t) — 
L\{t — z,t) : z £ (0,t)} are two independent Poisson processes. 

Define X'^{t) as the position at time t of a dual second class particle of 
the process Lx, that started in (0,0). Then we know that 

(5.2) x<y<X'^{t) =^ Lx{y,t)-Lx{x,t)<Lo{y,t)-Lo{x,t). 

This is due to the fact that if we leave out all the sources of Lx , the space- 
time paths do not change above the trajectory of X'^ (this is one of the key 
ideas in [4]; the reader might want to check this fact by looking at Figure 
2). This means that if we define the process L as a Hammersley process that 
uses the same a-points and sinks as Lx, but starts without any sources, we 
have that 

x<X;(t) =^ Lx{x,t) = L{x,t). 

Inequality (5.2) now follows from the fact that the set of particles of L is at 
all times a subset of the set of particles of the Hammersley process Lq, since 
this process has no sinks, whereas L does have sinks. 



Theorem 5.1. FixL>0. Then 



limsuppj sup {N{z)+At{z)}-At{0)>Lt^/^\=O{e^), e j 0. 



Proof. We will use the auxiliary Hammersley process constructed above. 
Define 

If X'^{t) > t, (5.2) tells us that, for < z < t, 

Lx{t, t) -Lx{t- z, t) < Lo(t, t) -Lo{t- z, t). 
Using (5.1) and defining 

Nx{z)=Lx{t,t)-Lx{t-z,t), 
we see that (remember that N and At are independent) 

p( sup {N{z) + At{z)}-At{0)>Lt^/A 

Ue[0,ei2/3] J 

(5.3) 

<p( sup {N{z)-Nx{z)}>Lt^/^]+F{X'x{t)<t}. 

UG[0,et2/3] J 
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The second term on the right-hand side of (5.3) can be bounded using 
(2.5), (3.5) and Corollary 4.4: 

F{X'^{t) <t}= F{X'{t/X) < Xt} 

= P{Z'(t/A)>t(l/A-A)} 

(5.4) <P{Z(t/A) >t(l/A- A)} 

= P{Z(t/(l - rt~^/^)) > ^2/3(2 - rt~^/'^)/{l - rr^/^)} 
<P{Z(t)>rt '/'}<r-3, 

for all r G [l,t^/^), applying Corollary 4.4 with argument i= t/{l — rt"^/^). 

The first term on the right-hand side of (5.3) concerns a hitting time for 
the difference of two independent Poisson processes. After rescaling, this can 
be written as 

P{3o<.<e : t-i/3{iV(zt2/3) _ N^(zt'/'^)} > L}. 

The process z i— > t~^^^{N{zt'^^^) — Nx{zt'^^^)} converges, as t ^ oo, to the 
drifting Brownian motion process 

(5.5) Wr{z)=^W{2z)+rz, z > 0, 

in the topology of uniform convergence on compacta, where W is stan- 
dard Brownian motion on R_|_. Hence, we get, by a standard application of 
Donsker's theorem, 



hm P{3o<.<, : t-'/'{N{zt^/') - Nx{zt^/')} > L} 

(5-6) 

sup Wr{z)>L\. 

GfO,el J 



2G[0,e; 

We now get, for r < L/e, 



I sup Wr(^)>^j<ip( sup W{z)>L-e 

Ue[0,el J U6[0,2e] 

sup W{2ez)/\/2e>{L-£r)/V2e\ 
6 [0,1] J 

\ sup W{z)>{L-er)/V2e 
Ue[o,il 



(5.7) 

= P<! £ 

z[0. 



(L-er)/V2£ 
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Taking r = L/{2e), we get 
'2 



vr J(L-sr)/V2^ 

g-{i/2)n2 2^/8ie-^V{i6.) 

, — dur^ = , ej,0, 

L/v^ V27r Lv27r 

using Mills' ratio approximation for the tail of a normal distribution in the 
last step. This means that, with this choice of r, our estimate for the second 
term on the right-hand side of (5.3) is dominant, so (5.4) now proves the 
theorem. □ 



Note that we could use the proof of the theorem to show that L can even 
go to at a certain speed when e — > 0, and still the considered probability 
would go to 0, uniformly in t. 

6. Lower bound for E.Z(f)_j-. We wish to bound the probability that 
Z{t) e [0,et2/3], In order to do this, we again introduce an independent 
auxiliary stationary Hammersley process Lai but now with intensity A > 1. 
In fact, we will choose 

Coupled to this process, we again consider the Hammersley process Lq with- 
out sources or sinks, but with the same a-points. This time, however, we 
will leave out the sinks of the stationary process; to be more precise, we have 
that 

(6.1) y>x>X),{t) =^ Lo{y,t)-Lo{x,t)<Lx{y,t)- Lx{x,t). 

The reason is that if we consider the stationary process Lx and leave out the 
sinks of this process, below the trajectory of Xx the space-time paths do not 
change. The inequality then follows from the fact that the set of particles of 
a process which, like Lq, has no sinks, but starts with sources and uses the 
same a-points, will at all times be a superset of the set of particles of the 
Hammersley process Lq. Compare this to the explanation of (5.2). 
We again define 

Nx{z) = Lx{t,t) - Lx{t - z,t) 
and will show the following result. 



Theorem 6.1. 



limlimsupP{0 < Z{t) < ei^/^} 
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Proof. Let 77 > 0. It is enough to find e > such that 
hmsupP{0 < Z{t) < et^/^} < Sry. 

t—*oo 

For any L, r > and A = 1 + rt~^^^, we have 

< < £^2/3} < p| sup (^N{z) + At{z))< sup (Af(z) + ^t(z))| 

U>£t2/3 ze[0,rt2/3] J 

(6.2) <p| sup {Niz)-{AtiO)-Atiz))}<Lt^/'^ 



+ p( sup {Niz)-{AtiO)-Atiz))}>Lt'/^\, 
Uefo.rf2/3i J 



since for any L G M, X < y imphes that either X < L or y > L, so 

F{X <Y)< F{{X < L} U {y > L}) < F{X <L) + P(L < y), 

which allows us to "optimize" over L. For the first term of (6.2), we want 
to use (6.1), which implies that, on the event {X\{t) <t — rt"^^^}, we know 
that, for ah z < rt^/^, 

Lo{t, t) -Lo{t- z, t) < Lx{t, t) -Lx{t- z, t). 

Therefore, 
P{0 < Z{t) < ef'^} 

4 sup {N{z) - Nx{z)} < Lt^/^l + F{Xx{t) > t - rt2/3} 



< 



(6.3) +P( sup {N{z)-{At{0)-At{z))}>Lt'/A 

UG[0,et2/3] J 

= p( sup {N{z)-Nx{z)}<Lt^^^]+F{X{t/X)>Xt-Xrf/^} 

Ue[ei2/3,ri2/3] J 

{N{z)-{At{0)-At{z))}>Lt^/''y 



sup 

2G[0,et2/3] 

Using (3.5) (note that At - Art^/s > t/X when r < ^t^/^), Z'{t) < Z{t) and 

—Z'{t) = Z{t), we get, for the second term in (6.3), 
F{X{t/X) >Xt- Art2/3} = w{Z{Xt - Xrt^'^) < (1/A - X)t + Art^/S} 

< P{Z'(At - Art2/3) < (1/A - X)t + Xrt^l^} 
= P{-Z'(At - Art2/3) > (A - l/A)t - Art2/3} 
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= F{Z{Xt - Xrt^/^) > (A - l/\)t - Art2/3} 

This means that we can start by choosing r sufficiently large to ensure that 
the second term is smaller than 77, since 



p|z(t(l-r2r2/3))>- 



2.1/3 



= p|z(t + 0{t^/^)) > rt2/3 + 0(ti/3)| ^ o(r-3), t ^ 00, 

where we use Corollary 4.4 in the last step. 

Now we turn to the first term in (6.3). This term is very similar to the 
term we found in the previous section. We have, as in the proof of Theorem 
5.1, that the process 

converges, as t ^ cxo, to a drifting Brownian motion process 
(6.4) Wr{z) =f W{2z) -rz, z> 0, 

in the topology of uniform convergence on compacta, where W is standard 
Brownian motion on (this time the drift is negative instead of positive) . 
Hence, we get, again using Donsker's theorem. 



(6.5) 



limp( sup {N{z) - Nx{z)} < Lt^/^\ 

UG[rf2/3,-rt2/3] J 

sup Wriz) <l\ 
£[e,r] ) 

sup Wr{z) <l\ + f\ sup Wr{z)<L, SUp Wr{z) > l\ 



e[0,r] J lze[e,r] ze[0,e] 

<P( sup Wr{z) <l\ + f\ sup Wr{z)>L 
UG[0,r] J Ue[0,e] 



Since 



limpj sup Wr{z) < L 1 = 0, 



we can choose L = L{rj) > sufficiently small to ensure 



j sup Wr{z)<L\ <'q/2. 
Ue[0,H J 
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It is also clear from the argument of the proof of Theorem 5.1 [see (5.7)] 
that we can next choose e > sufficiently small to ensure that 

p| sup Wr{z) > l1 < p( sup W{2z) >l\ < ?]/2, 

for this choice of L = L{r]) > 0. 

It is now seen from (6.5) that this bounds the first term of (6.3) from 
above by rj [remember that we have already fixed r > to bound the second 
term of (6.3)]. Finally, we can choose e > so small that the third term in 
(6.3) is smaller than rj, using Theorem 5.1. This completes the proof. □ 



Corollary 6.2. Let Z{t)+ and L{t,t) be defined as in (3.6). Then 

. ,EZ(t)+ ^. . „Var(L(t,t)) ^ 
hmmf — -^T-^ = hmmf ^, ' > 0. 

t^oc ^2/3 t^oo 2*2/3 



Proof. Using (3.6), we only have to prove the statement for EZ{t)^. 
Suppose t„ ^ CO such that 

2/3 ^ 



Then 



Since -Z'{t) = Z{t) and Z' {t) < Z{t), we have that F{Z{t) > 0} > 1/2, for 
all t>0. This would mean that, for all e > 0, 

liminf P{0 < Z(t„) < ei^/^} > i, 
which would contradict Theorem 6.1. □ 



Remark. We can again make a comparison with the results in [7]. Jo- 
hansson shows that the probability that all longest strictly NE paths from 
(0,0) to {t,t) stay within a strip around the diagonal of width t"' does not 
tend to 1, as t — > oo, for all 7 < 2/3, again using the analytical results in 
[2]. Our Theorem 6.1 shows that the probability that all longest weakly NE 
paths from (0,0) to {t,t) stay within a strip around the diagonal of width 
t"' tends to 0, as t — > 00, for all 7 < 2/3 (see also the Remark at the end of 
Section 4). 
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7. Tightness results. In the preceding sections it was shown, using the 
"hydrodynamical methods" of [4] that, for a stationary version of Hammer- 
sley's process, with intensity 1 for the Poisson point processes on the axes 
and in the plane, the variance of the length of a longest weakly NE path 
L{t,t) is of order t^/'^, in the sense that 

(7.1) < liminf Var(L(t,t)) < limsupi^^/^ Var(L(t, t)) < oo. 

This means, in particular, that, for any t > 0, the sequence 
n~^^^{L{nt,nt) -2nt}, n = l,..., 

is tight. 

As noted in [9], the distributional limit result for n~^^^{LQ{nx,nt) — 
2n\/xt} for Hammersley's process without sources and sinks in [2] can be 
translated into a limit of Yn/n^^^, where 

(7.2) Yn = Znt{[nx])-nxy{At), 

and Znt{[nx]) is the [nx]th particle at time nt, counting particles at time nt 
from the left. Theorem 3.2 in [9] gives a tightness result for a more general 
version of 5^, in the context of a version of Hammersley's process on the 
whole line, with a (possibly) random initial state. The result is that, under 
his conditions D and E, the sequence 

y„,/(n^/^logri), n = l,... 

is tight. He conjectures that, in fact, n^^^logn can be replaced by n^^^. 
The results, derived above, are a further indication that indeed n^/^logn 
might be replaced by n^/^, and that this can be derived by hydrodynamical 
methods. 

For the stationary version of Hammersley's process, with intensities 1 of 
the Poisson processes in the plane and on the axes, we can define Znti[2nt]) 
as the location of the [2nt]th source at time nt, where we count the sources 
from left to right, starting with the first source to the right of zero. Note 
that at time zero the particles are just the sources. The particles, escaping 
through a sink, are given location zero at times larger than or equal to the 
time of escape. 

With this definition, our results give tightness of the sequences 

(7.3) n-^/'^{znt{[2nt])-nt}, n = l,2,..., 

for each t > 0. This can be seen in the following way. We have, for Af > 
and t > 0, the "switch relation" 

(7.4) n~'^^^{znt{[2nt]) -nt}> M L{nt + Mn^^^ ,nt) < [2nt]. 
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Theorem 2.1 yields 

Yai{L{nt + Mn^/^,nt)) = -nt - Mn^/^ + nt + 2E(nt + Mn^/^ - Xi{nt))^ 

= -Mn^/^ + 2E(nt + Mn^'^ - Xi{nt))^ 
< Mv}/'^ + 2E(nt - 

where we use {Y + Z)+ < y+ + Z+ in the last step. Theorem 2.1, applied in 
the opposite direction, yields 

2R{nt - Xi{nt))^=\ai{L{nt,nt)). 

Hence, we get, by (7.1) and Chebyshev's inequality, 

¥{n-^/^{zni{[2nt\) - nt] > M] 

= W{L{nt + Mn^/^, nt) - 2nt - Mn^/^ < [2nt] - 2nt - Mn^/^} 

MnV3+Var(Li(nt,nt)) ^ Mn^/^ + 0{{ntf/^) _ ^, ^ 
- {Mn^/^ + 2nt-[2nt]Y ~ NWT^ ~ ^ 

We similarly get 

f{n~'^/^{znt{[2nt\) - nt} < -M} = 0(M-^), 
using that, if nt — Mn~^/^ > 0, 

(7.5) n^^/^{znt{[2nt]) - nt] < -M L{nt - Mn^'^ ,nt) > [2nt], 

which proves the tightness of the sequences (7.3). 

Although the tightness of the sequence (l^/n^/^) for the Hammersley 
process without sources or sinks, as defined in (7.2), is known from the 
results of [2], it is of some interest to derive this from the results of the 
preceding sections. The tightness will follow from 

(7.6) n^^/^{Lo{nx,nt) -2y/xi} = Op{l), n ^ oo, 

for all x,t>0, where LQ{nx,xt) is a strictly NE path from (0,0) to (nx,xt), 
and where the intensity of the Poisson process in the first quadrant is equal 
to 1. 

We again have a "switch relation" similar to (7.4): 

n~'^/^{znti[2nt]) -nt} > M Lo{nt + Mn^^^ ,nt) < [2nt]. 

From [1], we know that 

Lo(nx,nt) = LQ{n\^,ny/xi), 
so if we show that, for each t>0, 

(7.7) n~'^/^{Lo{nt,nt) -2nt} = Op{l), n^oo, 
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we get, for each e > and i > 0, 
P{n~^/^{z„t([2nt]) - nt} > M} 

= F{Lo{nt + Mn^^^,nt) - 2nt - Mn^^^ < [2nt] - 2nt - Mn^/^} 

= F{Loint^Jl + Mn-'^/^/t, nt^J I + Mn-"^/'^ /t) 
- 2nty^l + Mn-2/3/t + 0{M'^n-^/^) 

< [2nt] - 2nt - Mn^''^} 

<e, 

for sufficiently large M = M(e) > and all n > no(A/, e). Relation (7.7) 
similarly implies 

¥{n-^/^{znt{[2nt]) - nt} < -M} < e, 

for sufficiently large M = M{e) > and all n > no(-/Vf, e), using 

n~'^/^{znt{[2nt]) - nt} < -M Lo{nt - Mn^/^,nt) > [2nt]. 

In order to prove (7.7), it is sufficient to show 

(7.8) t^^/3{Lo(t,t) - 2t} = Op(l), t-^oo. 

Now first note that the length LQ{t,t) of a longest strictly NE path from 
(0,0) to {t,t) is the same as At{0) in the proof of Theorem 5.1. Let, for 
A = 1 — rt~^^^, L\ be defined as in the proof of Theorem 5.1. By (5.3), we 
have 

p( sup {N{z)+Atiz)}-At{0)>Lt^/^] 

UG[0,/^t2/3] J 

(7.9) 

<p( sup {Niz)-Nx{z)}>Lt^/^]+F{X'^it)<t}. 
Ue[o,i^t2/3] J 

We first deal with the second term on the right-hand side of (7.9). By (5.4), 
we have 

(7.10) F{X',{t)<t} = 0{r-^), 

uniformly for r E [1, ^t^^'^]. To deal with the first term on the right-hand side 
of (7.9), we first note that 

z ^ Mr{z) =^ t~^/^{N{zt^/^) - Nxizt^/^)} -rz, z> 0, 
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is a zero-mean martingale. So we get 

p( sup {N{z)-Nxiz)}>Lt^/^\=¥{ sup {Mriz)+rz}>L 

<p| sup Mr{z)>L-rK\. 
Taking r = L/[2K) and using Doob's submartingale inequality, we get 



p( sup {N{z)-Nx{z)}>Lt^'^\<A sup Mr{z)>L/2 

<p( sup Mr{zf>L'^/A] 



We also have, by Corollary 4.4, 

P( sup {Niz)+At{z)}^ sup {iV(z)+^t(z)}j 

< 2P{Z(t) > i^t2/3| < ;^/^3_ 

So, taking K = l?!^"^ [note that this means that r = L/{2K) = \L^I^'^ < 
for L < t^/S], we obtain 

p( sup {N{z) + Atiz)}-At{0)>Lt'/'^] 
Ue[-t,t] J 

= P{Li(t,t)-^i(0)>Lti/3}<L-5/4^ 

for all L<t4/5. 

If L > t^/^ , we first note that 

F{Li{t,t) - AtiO) > Lt^/^} <F{Li{t,t) > Lt^/^} 

< 2F{Pt > ^Lt^^^} < 2F{Pt > ^L^^^t}, 

where Pt is a Poisson variable with expectation t. Let [x] denote the largest 
integer < x and let a = ^L^^^. Then, using the Lagrange remainder term in 
an expansion of e*, we get, for a 9 £ (0, 1), 

.[at] -{l-e)t Jat] 

P{P,>lLV.„<P{P.>M} = ^^<pp. 
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Stirling's formula for the gamma function r(x) yields that, uniformly in 
t> 1, 

j.at 1 

^g-a(loga-l)t 



T{at+1) 

This implies that P{Pt > ^Lt^^^} tends to zero faster than any negative 
power of L, if L > t^/^, uniformly in all large t and hence, we can conclude 
that 

P{Li(t,t) - MO) > = 0(L-5/4), 

for all L> 1 , implying 

0<2t-EAt{0)=lE{Li{t,t)-AtiO)} 

(7.11) 
Thus, 

E\At{0) - 2t\ = E|Lo(t, t) - 2t\ = 0{t^/^). 

This proves (7.8) and, as noted above, (7.6) now also follows. This, in turn, 
proves tightness of the sequence {Yn/v}/^), for y„ defined by (7.2) for Ham- 
mersley's process, starting with the empty configuration on the axes. 
Notice that we also proved at the same time 

ELQ{nx,nt) = ELQ{n^/xi,nVxi) = EA^^^{Ql) = 2nVxt + 0{n^^^), 

for all x,t>0; see (7.11). 
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